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Abstract 

We work out a superstatistical description of high-energy scat- 
tering processes that takes into account temperature fluctuations in 
small volume elements. For F-distributed fluctuations of the inverse 
temperature one effectively obtains formulas similar to those used in 
nonextensive statistical mechanics, whereas for other temperature dis- 
tributions more general superstatistical models arise. We consider two 
main examples: Scattering processes of cosmic ray particles and e^e~ 
annihilation processes. In both cases one obtains excellent fits of ex- 
perimentally measured energy spectra and cross sections. 



1 Introduction 



Super statistical techniques have been recently successfully applied to a large 
variety of complex systems, for example hydrodynamic turbulence [H 121 [3l H] , 
defect turbulence [5] , share price dynamics [HI [7j , random matrix theory [8], [9] , 
random networks ^Q\, wind velocity fluctuations [HI [12], hydro-climatic fluc- 
tuations [13], the statistics of train departure delays [H] and models of the 
metastatic cascade in cancerous systems [15]. The basic idea underlying 
this approach is that there is an intensive parameter, for example the in- 
verse temperature /3 or the energy dissipation in turbulent systems, that 
exhibits fluctuations on a large time scale (large as compared to internal re- 
laxation times of the system under consideration). As a consequence, one 
can model these types of complex systems by a kind of superposition of ordi- 
nary statistical mechanics with varying temperature parameters, in short a 
superstatistics [l6l[I7l[T8l[T9l[20l[211[22l[23]. The stationary distributions of 
superstatistical systems deviate from ordinary Boltzmann-Gibbs statistical 
mechanics and can exhibit asymptotic power laws, stretched exponentials, or 
other functional forms in the energy E [18j. 

In this paper we work out potential applications of this concept in high 
energy physics. In scattering processes at high energies, the effective inter- 
action volume, as well as the number of particles involved, can be rather 
small, and hence temperature fluctuations can play a very important role 
[2H [25l [26] . These can either be temperature fluctuations from scattering 
event to scattering event, i.e. temporal fluctuations, or they may be related 
a nonequilibrium situtations where different spatial regions have different 
temperature, i.e. spatial fluctuations. 

Our starting point for a suitable thermodynamic model will be Hagedorn's 
theory [271 128], which yields a statistical description of a selfsimilar 'fireball' 
of particles produced in scattering events. Hagedorn's theory models the 
hadronization cascade from a statistical mechanics point of view (of course 
QCD was not known at the time when he wrote the seminal paper [2Z])- His 
theory is regularly in use to describe heavy ion collisions [23 [SHI [SI]- The 
basic assumption is that in the scattering region the density of states grows 
so rapidly that the effective temperature cannot exceed a certain maximum 
temperature, the Hagedorn temperature Th [32]. The value is approximately 
Tfj ^ 180 MeV and it describes the confinement phase transition. The 
Hagedorn phase transition is also of fundamental interest in string theories 

[231 EH [32]. 
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In this paper we will consider a superstatistical extension of the Hage- 
dorn theory, starting from a g-generalized version previously introduced in 
[36]. The predictions of this generalized statistical mechanics model are in 
excellent agreement with measured experimental data. We will illustrate 
this for two main examples: i) observed energy spectra of cosmic rays [21] ii) 
experimentally measured cross sections in e^e~ annihilation 

The Hagedorn theory of scattering processes is known to give correct 
predictions of cross sections and energy spectra for center of mass energies 
EcMS < ^OGeV, whereas for larger energies there is experimental evidence 
from various collision experiments that power-law behaviour of differential 
cross sections sets in. This power-law is not contained in the original Hage- 
dorn theory but can be formally obtained if one extends the original Hagedorn 
theory to a superstatistical one, along the lines sketched in this paper. In 
fact, if fluctuations are taken into account then power laws can arise in quite 
a natural way in various ways, and often lead to Tsallis type of generalized 
distribution functions[37l [38l [39l [H |16]. While in this paper we concentrate 
on cosmic ray statistics and transverse momentum spectra in e~^e~ annihila- 
tion, it is quite clear that related superstatistical techniques can be applied 
to other scattering data as well, for example heavy ion collisions [311 SO] , PP 
collisions [H] and pp collisions [26] . 

2 Superstatistical modelling of temperature 
fluctuations in scattering experiments 

In the superstatistics approach [16] one assumes that locally the system un- 
der consideration reaches local equilibrium but on a large spatio-temporal 
scale there are temperature fluctuations, described by a probability density 
/(/?) which describes the probability to observe a certain inverse temperature 
/3 in a given spatial area. One thus gets a kind of mixing (or superposition) of 
many equilibrium distributions which effectively decribe the driven nonequi- 
librium system with a stationary state (see [12] for a recent review discussing 
various applications). 

For high energy scattering processes, it is clear that the larger the center 
of mass energy Ecms of the collision process is, the smaller is the volume 
probed, due to the uncertainty relation \Ecms " ~ 0{h). This means, 
the effective interaction volume where a thermodynamic description of the 
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collision process makes sense will become smaller and smaller with increasing 
EcMS- However, a smaller volume means larger temperature fluctuations. 
This is in particular relevant if we repeat our scattering experiment several 
times or if we have different scattering events in different spatial regions 
which all contribute to our data. It thus makes sense to consider at large 
energies Ecus a superstatistical description of scattering events which takes 
into account local temperature fluctuations. 

Assume that locally some value of the fluctuating inverse temperature j3 
is given. We then expect the momentum of a randomly picked particle in 
this region to be distributed according to the relativistic Maxwell-Boltzmann 
distribution 

P{E\P) = zWf 

Here p{E\(3) denotes the conditional probability to observe a particle with 
energy E = ^J-p^c^ + m^c^, given some value of /?. For highly relativistic 
particles we can neglect the rest mass m so that E = c\p\, where p is the 
momentum. The normalization constant is given by 

Z{P) = E\-^^dE = ^^. (2) 

Now let us take into account local temperature fluctuations in the small 
interaction volumes where scattered particles are produced. We have to 
consider some suitable probability density f{l3) of the inverse temperature in 
the various interaction volumes. In a long series of experiments we will then 
observe the marginal distribution obtained by integrating over all /3 

/>oo 

p{E) = / p{E\P)f{P)dp. (3) 
Jo 

These types of distributions are generally studied in superstatistical models, 
and are known to exhibit fat tails whose asymptotic decay with E depends 
on how f{j3) behaves for /? — > (more details on this in [IHj). Of course, 
the validity of the superstatistical modelling approach requires sufficient time 
scale separation so that the system relaxes to local equilibrium fast enough 
before the next temperature fluctuation takes place [31 [71 [201 123] . 

It has been argued in previous work [3] that there are basically three 
different distributions /(/?) that are relevant for large classes of complex sys- 
tems: These are the x^"distribution, the inverse distribution, and the 
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lognormal distribution. Lognormal superstatistics is often observed in hy- 
drodynamic turbulence, due to the multiplicative nature of the Richardson 
cascade [H H] . Inverse superstatistics plays an important role in random 
matrix theory [HI [9], as well as in medical statistics [15]. For high energy 
physics, superstatistics is most relevant, although other superstatistics 
may play a role as well (for example, one could think about a superstatistics 
generating Kaniadakis statistics [15]). 

For superstatistics (or equivalently F superstatistics), the distribution 
/(/?) is given by the x^-distribution of degree n, i.e. 

The x^-distribution is a typical distribution that naturally arises in many cir- 
cumstances, for example if n independent Gaussian random variables Xj, z = 
1, . . . , n with average and the same variance are squared and added. If we 
write 

fi--=Y.n (5) 

i=l 

then /? has the probability density function (j4]). The average of the fluctuating 
/3 is given by 

{/3) = n{Xf) = / (3fi/3)d(3 = /3o (6) 



and the variance by 

m-p', = -Pi (7) 

n 

The integral ([3]) with /(/5) given by (jl]) and p{E\i3) given by ([1]) is easily 
evaluated and one obtains 

p{E) (8) 

(1 + %- 1)E)^ 



where 



and 
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Note that the partition function Z{P) entering into eq. is /3-dependent. 
Different /^-dependencies of the partition function Z{f]) lead to different an- 
swers if the integration over (3 is performed. In particular, the precise relation 
between q and n depends on this. This was for the first time correctly worked 
out in [1]. 

The distribution ([H]) is a g-generalized relativistic Maxwell-Boltzmann dis- 
tribution in the formalism of nonextensive statistical mechanics These 
kind of distributions can be directly obtained by maximizing the g-entropies 



and multiplying with the available phase space volume. The pi are the prob- 
abilities of the microstates i. The g-entropies contain the Shannon entropy 
5*1 = —k^^pilnpi underlying ordinary statistical mechanics as a special 
case for q = 1. Whereas g = 1 corresponds to the usual canonical en- 
semble with constant temperature, the Tsallis-canonical ensemble obtained 
for g > 1 is capable of describing temperature fluctuations, assuming that 
P is x^-distributed. We thus have a plausible physical mechanism why a 
Tsallis-like statistical description makes sense if a suitable intensive param- 
eter fluctuates [U [38]. Tsallis statistics with g = 3 (respectively g = — 1 if 
the escort formalism HSj is used) also plays an important role for chaot- 
ically quantized scalar fields [l6l HT], dark energy models jH] and moduli 
field dynamics ^9] . 

3 Superstatistical partition functions 

When experimental scattering data are collected, one often looks at momen- 
tum distributions of a particular particle that are produced by repeating the 
same experiment many times. In each scattering event, the effective temper- 
ature (given by the heat bath of surrounding particles) will fluctuate from 
event to event. Since the scattering experiment is repeated many times in 
an independent way, and since we concentrate on the statistics of just one 
particle rather then many-particle states, it makes sense to directly integrate 
out the temperature fluctuations and to consider effective 1-particle super- 
statistical Boltzmann factors. 

Let us consider particles of different types and label the particle types by 
an index j. Each particle can be in a certain momentum state labelled by 



m 




(11) 
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the index i. The energy associated with this state is 



^Jpf + ^l (12) 



using units where c = 1. We may now define an effective 1-particle Boltzmann 
factor bjQ 

/■oo 

X,, := / rf/?/(/3)e-^^-. (13) 
Jo 

For example, for Tsalhs statistics /(/?) is a distribution and one has 

Xi,- = (l + (g-l)6ei,)-^, (14) 

where is proportional to the average temperature, and g — 1 is a measure 
of the strength of inverse temperature fiuctuations. The effective Boltzmann 
factor Xij approaches the ordinary Boltzmann factor e"**^'^ for g — 1. 

The simplest way to introduce a generalized grand canonical partition 
function (which, as said before, regards all temperature fluctuations to be 
independent) is as follows: 

Here denotes the number of particles of type j in momentum state i. The 
sum stands for a summation over all possible particle numbers. For 

bosons one has Uij = 0, 1, 2, ... , oo, whereas for fermions one has i>ij = 0, 1. 
It follows that for bosons 



^^x^j-' = (bosons) (16) 

1 Xij 

whereas for fermions 



x^j-' = 1 + rCjj (fermions) (17) 

Hence the partition function can be written as 

Z = l[-^l[(l + Xiy), (18) 

-L Xi^ 



^To describe 2-particle states with the same temperature fluctuations surrounding both 
particles one has to proceed in a slightly different way, see [1]. 
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where the prime labels fermionic particles. For the logarithm of the partition 
function we obtain 

logZ = - J] log(l - Xij) + J2 log(l + ^i'j')- (19) 

ij i'j' 

One may actually proceed to continuous variables by replacing 
(Vq: volume of the interaction region) and 

/•oo 

- / P{m)[...]dm, (21) 



where p(m) is the mass spectrum. 

Let us now formally calculate the average occupation number of a particle 
of species j in the momentum state i. We obtain 

^'^^=^^^BZ-.^^^Z = -^. (22) 
For example, for the case of TsaUis statistics one obtains 

= ^ 1 (23) 

(l + (g-l)6e,,)5^±l 

where the — sign is for bosons and the + sign for fermions. 

In order to single out a particular particle of mass mo, one can formally 
work with the mass spectrum p{m) — 5{m — mo). To obtain the probability 
to observe a particle of mass ttlq in a certain momentum state, we have to 
multiply the average occupation number with the available volume in mo- 
mentum space. An infinitesimal volume in momentum space can be written 
as 

dpxdpydpz — dpLPT sin OdprdO (24) 

where pt = a/p^ + pI is the transverse monentum and Px = Pl is the longi- 
tutinal one. By integrating over all 9 and pl one finally arrives at a proba- 
bility density w{pt) of transverse momenta given by 

/■cx> Y 

w{pt) = const ■ Pt / dpi j . (25) 

^0 (l + (g- l)6v'p2 +m^)^ ±1 



Since the Hagedorn temperature is rather small (of the order of the vr mass) , 
and is of the order of the Hagedorn temperature, under normal circum- 
stances one has b^/p"^ + pt + "^o 1' hence the ±1 can be neglected if 
q is close to 1. One thus obtains for both fermions and bosons the statistics 



w{pt) ~ const ■ pt j dpi yl + {q — l)h^ p^ + p| + mgj (26) 

which, if our model assumptions are satisfied, should determine the pr de- 
pendence of experimentally measured particle spectra. The differential cross 
section a'^dajdpT is expected to be proportional to w{pt)- 

In nonextensive statistical mechanics, it is sometimes of advantage to 
consider normalized g-expectation values [35], thus implementing a formal- 
ism that is based on so-called escort distributions [44j. If the escort formalism 
is used then the power the above formula is replaced by ^^y, a simple 

reparametrization. For more general classes of superstatistics described by 
general /(/?) the g-exponential is replaced by /(/3) exp{—(3^^/p'^ + pt + ml}d(3. 

4 Comparison with measured cosmic ray en- 
ergy spectra 

Let us now look at measured data. Fig. 1 shows the experimentally measured 
energy spectrum of primary cosmic ray particles as observed on the earth 
[50l ISTl [52l [53l EH [55]. The measured spectrum is very well fitted over many 
decades of different energies by the distribution ([8]) [24J or other generalized 
nonextensive distributions [56]. The best fit is obtained if the entropic index 
q is chosen as 

g= 1.215 (27) 
and if the effective temperature parameter is given by 

kTo = b-^ = 107 MeV. (28) 

Let us now predict a plausible value of q using the superstatistical ap- 
proach. The variables Xj in eq. describe the independent degrees of 
freedom contributing to the fluctuating inverse temperature (3. At very large 
center of mass energies Ecms oo, the interaction region is very small, 
and all relevant degrees of freedom are basically represented by the 3 spatial 
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Figure 1: Measured energy spectrum of primary cosmic rays (in units of 
m~^s~^sr~^GeV~^) [23]. The solid line is the formula ([8]) with q = 1.215, 
b^^ = kTo = 107 MeV and C = 5 • 10~^^ in the above units. The dashed line 
is eq. (EI) with q = 11/9, kTo = 107 MeV and C smaller by a factor 1/50. 
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dimensions into which heat can flow. We may physically interpret Xf as the 
heat loss in the spatial i-direction, i = x,y, z, during the collision process 
that generates the primary cosmic ray particle. The more heat is lost, the 
smaller is the local kT, i.e. the larger is the local (3 given by ([5]). The 3 
spatial degrees of freedom yield n = 3 as the smallest possible value of n or, 
according to 

g = y = 1.222. (29) 

For cosmic rays Eqms is very large, hence we expect a g- value that is close to 
this asymptotic value. The fit in Fig. 1 in fact uses q = 1.215, which agrees 
with the predicted value in eq. (jSHD to about 3 digits (similar g-values were 
also obtained in [5B]). 

For smaller center of mass energies, according to \Ecms ■ ~ 0{h), the 
interaction region will be bigger and more effective degrees of freedom within 
this bigger interaction region will contribute to the fluctuating temperature. 
Hence we expect that for smaller Ecms n will be larger than 3, or q will be 
smaller than 11/9. 

Finally, for Ecms 0, q ^ 1 and ordinary statistical mechanics is recov- 
ered. In this classical limit case, the relevant interaction region where a 
thermodynamic description makes sense becomes very large, and within this 
large region a large number of independent degrees of freedom n contributes 
to the fluctuating temperature, represented by many different particles. Ac- 
cording to eq. ([9]), the limit n — > cxd is equivalent to g — > 1 and in this limit 
the x^-distribution degenerates to a delta function 5(/5 — Pq). 

It is reasonable to assume [51] that the 'knee' at ~ 10^^ eV is due to the 
fact that one has reached the maximum energy scale to which typical galactic 
accelerators can accelerate. This then implies a rapid fall in the number of 
observed events with a higher energy, i.e. a steeper slope in Fig. 1 between 
about 10^^ and 10^^ eV. The 'ankle' at E fa 10^^ eV may then be due to 
the fact that a higher energy population of cosmic ray particles takes over 
from a lower energy population. This higher energy population may have 
a different origin (for example, extragalactic origin). The new population 
has a significantly smaller flux rate but can reach much larger energies. As 
a matter of fact the cosmic accelerators underlying the production process 
of this new species of cosmic rays must have a much larger center of mass 
energy Ecms than the ankle energy ~ 10^^ eV, so q should be given by its 
asymptotic value 11/9, whereas the effective temperature Tq should be the 
same as before. The dashed line in Fig. 1 corresponds to our formula with 
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E [oVJ 

Figure 2: Measured cosmic ray energy spectrum ■ dN/dE at largest en- 
ergies (data from [51], [521 [531 El])- The straight line is a power law with 
exponent a = 5/2 (corresponding to q = 11/9). 

q = 11/9, kT = 107 MeV and a flux rate that is smaller by a factor 1/50 as 
compared to the high-flux generation of cosmic rays. This is consistent with 
the data. 

Formula ([H]) predicts asymptotic power-law behavior of the measured en- 
ergy spectrum. For large E one has p{E) ~ E'"^ where the index a is given 
by 

a = — 2. (30) 

q-1 

q = 1.215 implies a = 2.65 (for moderately large energies), whereas for 
largest energies the asymptotic value q = 11/9 implies a = 5/2. As shown 
in Fig. 2, the largest-energy events are compatible with such an asymptotic 
power law exponent. 
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5 Comparison with experimentally measured 
differential cross sections in e^e~ annihila- 
tion 

Let us now proceed to our second example, differential cross sections for 
transverse momenta in e~^e~ annihilation experiments [361 EZ]- If one uses 
the escort formalism |^ HS] then formula fl26l) implies 



1 da 



/oo 2_ 
dx (^1 + (g - 1)^x2 + u2 + m2j (31) 



a dpT 

Here x = Pl/Tq, u = Pt/Tq and mp := uiq/Tq are the longitudinal momen- 
tum, transverse momentum and mass in units of a temperature parameter 
To that is of the same order of magnitude as the Hagedorn temperature. C 
is a suitable constant related to multiplicity. 

Generally, one knows that the interaction energy in a nonextensive system 
increases with increasing entropic index q. Since this energy must be taken 
from somewhere, it is most natural to assume that it is taken from the heat 
bath. So the average temperature parameter Tq should slighly decrease with 
increasing q. In [36] a linear dependence of Tq on q was postulated, of the 
simple form 

To = (l - I) Th (32) 

where Th = 180 MeV is the Hagedorn temperature. 

We also have to clarify the dependence of q on Ecus- Clearly, for 
EcMs oo our argument in section 2 and 4 suggests the value q = 11/9, 
whereas for Eqms ~^ one has g = 1, i.e the ordinary Hagedorn theory 
without fluctuations. In [36], a smooth interpolation between these values 
was suggested, of the form 

(i(^CMs) = 9 ^ ^-EcMs/Eo (33) 

where Eq ^ 45.6 GeV is about half of the Zq mass. 

Finally, one also needs to know the multiplicity (the average number 
of produced charged particles) in order to estimate the cross section. The 
multiplicity M as a function of Ecus has been independently measured in 
many experiments [5H]. A good fit of the experimental data in the relevant 



12 




0.1 1 10 



Figure 3: Differential cross section as a function of the transverse momentum 
Pt for various center of mass energies E. The data correspond to measure- 
ments of the TASSO {E < 44 GeV) and DELPHI {E > 91 GeV) collabora- 
tion. The solid lines are given by the analytic formula ( l35|) . 
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energy region is the formula 



M = (34) 

The final formula derived for the cross section in involves a further ap- 
proximation step to perform the integration over x and is then finally given 

by 

= ^Mp{u). (35) 
a dpT To 

where p{u) is the normalized probability density 

p{u) = ^ (1 + (g - 1)m)"^+^ (36) 
with normalization constant 

Z, = (,-irV^BQ,_^_3)^ (37) 
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Formula fl5^ with p{u) given by eq. fl5Bl) . q{EcMs) given by eq. 
To(g) given by eq. (l32|) and multiplicity M{Ecms) given by eq. (IMl) turns 
out to very well reproduce the experimental results of measured transversal 
cross sections for all energies Ecus- This is illustrated in Fig. 3 which shows 
the measured differential cross section versus pt- The solid lines are given 
by formula for the various center of mass energies Ecus- One obtains 
excellent agreement with the measured data of the TASSO and DELPHI 
collaborations [591 EO]- Remarkably, for the largest energies the best fitting 
parameter q is again given by g ~ 11/9, in agreement with the data for the 
cosmic rays. Moreover, at the largest energies, according to eq. fl52]) Tq is 
given by 107 MeV , again in agreement with what was used for the cosmic 
ray data in Fig. 1. Our superstatistical argument in section 2 and 4, relating 
q to temperature fiuctuatuons and predicting from this g 11/9 for largest 
energies Ecms, is indeed generally applicable. 



6 Conclusion 

In this paper we have dealt with a superstatistical generalization of the Hage- 
dorn theory which is generally applicable to describe the statistics of scattered 
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particles produced in high-energy colhsions. Our approach is based on taking 
into account temperature fluctuations in small interaction volumes. For x^- 
distributed inverse temperature one effectively ends up with formulas that 
are similar to those used in non-extensive versions of statistical mechanics. 

At large energies the x^"Superstatistical approach implies that energy 
spectra of particles and cross sections decay with a power law. This power 
law is indeed observed for various experimental data. We obtained formulas 
that are in very good agreement with experimentally measured data for cos- 
mic rays and e^e~ annihilation. In particular, the superstatistical approach 
allowed us to give a concrete prediction for q at largest center of mass en- 
ergies, namely q ~ 11/9. Generally, it appears that high-energy scattering 
data do not only yield valuable information on elementary particle physics, 
but they may also be regarded as test grounds to further develop generahzed 
versions of statistical mechanics. 
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